arXiv: 1507.0238lv3 [cond-mat.dis-nn] 11 Feb 2016 


Effect of Disorder in a Three-Dimensional Layered Chern Insulator 

Shang Liu,^ Tomi Ohtsuki,^ and Ryuichi Shindou^’^’0 

^School of Physics, Peking University, Beijing 100871, China 
^Department of Physics, Sophia University, Chiyoda-ku, Tokyo, 102-8554, Japan 
^International Center for Quantum Materials, Beijing, 100871, China 
^Collaborative Innovation Center of Quantum Matter, Beijing, 100871, China 

We studied effects of disorder in a three dimensional layered Chern insulator. By calculating the 
localization length and density of states numerically, we found two distict types of metallic phases 
between Anderson insulator and Chern insulator; one is diffusive metallic (DM) phase and the other 
is renormalized Weyl semimetal (WSM) phase. We show that longitudinal conductivity at the zero 
energy state remains finite in the renormalizd WSM phase as well as in the DM phase, while goes 
to zero at a semimetal-metal quantum phase transition point between these two. Based on the 
Einstein relation combined with the self-consistent Born analysis, we give a conductivity scaling 
near the quantum transition point. 
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Introduction. — During the last decade, huge research 
effort has been made on topological phases in condensed 
matter physics mu. Among many others, the criti¬ 
cal nature of the quantum phase transition between the 
topological phase and nontopological phase is one of the 
most fundamental research issues in the studies of topo¬ 
logical phases. A classic example is the localization prob¬ 
lem in the two-dimensional (2D) integer quantum Hall 
(QH) phase 0 (or so-called ‘Chern insulator’ 0-0), 
where a bulk delocalized state is universally observed 
between topologically distinct integer QH phases 
The existence of the bulk delocalized state comes from 
the stability of topological edge modes in respective in¬ 
teger QH phases 0, while its ‘zero-measure’ property 
is compatible with the scaling theory of the Anderson lo¬ 
calization ini- In fact, weak interlayer couplings among 
2D QH layers make the QH quantum critical point into 
a finite width of 3D metallic phase mm. 


In theoretical studies, preceding efforts demonstrated 
that multiple-layered Chern insulators (Cl) result in a 
novel semimetal phase, dubbed as ‘Weyl semimetal’ [ 3 - 


17j , where a pair of band touching points appear in the 


3D momentum space. The pair can be regarded as topo¬ 
logical defects having opposite magnetic charges, so to 
speak, ‘magnetic monopole’ and ‘antimonopole’. On in¬ 
creasing interlayer couplings, the monopole in relative to 
the antimonopole winds up the 3D Brillouin zone at least 
once, only to annihilate with the antimonopole again. 
This pair creation and annihilation process changes a 
global topological feature encoded in electronic Bloch 
wavefunctions, which ‘allows’ the system to enter into an¬ 
other kind of topological phase or non-topological gapped 
phase. When comparing these Weyl semimetal (WSM) 
physics with localization problem in integer QH physics, 
one may naturally ask ‘how the 2D QH quantum critical 
point is connected with 3D WSM phase ?’ or ‘Can we 
obtain a new insight to the localization physics in 2D QH 
system from the viewpoint of the 3D WSM physics ?’ 



FIG. 1. (color online) Phase diagram of layered Chern in¬ 
sulator subtended by disorder strength W and the interlayer 
coupling strength \f3\. WSM, DM, Cl, AI and 2dQHC stand 
for renomarlized Weyl semimetal, diffusive metal, Chern in¬ 
sulator, ordinary Anderson insulator and 2D QH quantum 
critical point, respectively, (inset) Phase diagram obtained 
from the self-consistent Born analysis. The region-II WSM 
phase (WSM-H) has two pairs of the Weyl points, while the 
region-III WSM phase (WSM-III) has three pairs of the Weyl 
points la. 


In this Letter, we will answer these questions, by study¬ 
ing numerically disorder effects in a layered Chern in¬ 
sulator. A computation of the localization length and 
two-terminal conductance by the transfer matrix (TM) 
method [isl and density of states by the kernel polyno¬ 
mial expansion (KPE) method [l9| reveals two distinct 
kinds of metallic phases, which intervene between con¬ 
ventional Anderson insulator (AI) and layered Chern in- 
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sulator (Cl) in a phase diagram subtended by the layer 
coupling and the disorder strength; one is (i) diffusive 
metallic (DM) phase with its zero-energy state, i.e. an 
electronic state at the band touching point (E* = 0), hav¬ 
ing a finite life time and the other is (ii) renormalized 
WSM phase whose zero-energy state has an infinite life 
time while the velocity and location in the k space of gap¬ 
less Weyl fermions being strongly renormalized by disor¬ 
der. A scaling analysis of the density of states clarifies a 
critical nature of the phase transition line between these 
two metallic phases. A finite-size-scaling of the conduc¬ 
tance also shows that finite conductivities in DM phase 
and in WSM phase go to zero toward the quantum crit¬ 
ical line. 

Model - We study a spinless two-orbital tight-binding 
model on a cubic lattice, which comprises of s-orbital and 
p = Px-\- ipy orbital; @ 


^ = E ([^* + ^'^(®)]4,sC£c, 5 + [ep + Vp(a;)]ctpCa,,p) 

X 

+ E(“ E 




X iJ,=x,y 


+ tsp 
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-itsp^ic 


o) ^x,s 


{'^s^x-\-ez,s^x,s + '^p^x-\-ez,p^^^P. 


+' J 


,)+H.c.), (1) 


where e^, Cp and Vs{x)^ Vp{x) denote atomic energies for 
the s, p orbital and disorder potential for the s, p or¬ 
bital, respectively. Both Vs{x) and Vp{x) are uniformly 
distributed within [—W/2,W/2] with identical proba¬ 
bility distribution. E, tp and tgp are intralayer trans¬ 
fer integrals between neighboring s orbitals, p orbitals 
and that between s and p orbital, respectively, while 
t'g and tp are interlayer transfer integrals. Without in¬ 
terlayer coupling, the model reduces to a 2D Cl, given 
that the so-called band inversion condition is satisfied; 
0 < \es — ep\ < 4:{ts tp). In the rest of this Letter, 
we take = —2{ts + tp), while t'^ = —tp > 0 

and ts = tv> 0, tsp = 4ts/3 with Ms being the en¬ 
ergy unit. [20| This parameter set realizes a Cl with a 

large band gap in the 2D limit {/3 = ^ = 0) for 

the half filling case. When 0 < |yd| < 1/2, the system 
is fully gapped, which belongs to the Cl phase in the 
2D limit. When |/3| > 1/2, the system enters into 3D 
WSM phase, where three pairs of the monopoles and an- 
timonopoles appear ai k = (0, tt, 7r±/co), (tt, 0,7r±/co) and 
(tt, TT, =F^i) respectively. The Weyl points at (0, tt, tt —/cq), 
(tt, 0, TT — /co) and (tt, tt, —/ci) are the magnetic monopoles 
with positive magnetic charge, while the others are the 
antimonopoles with negative charge. At \/3\ = three 
pairs of monopoles and antimonopoles annihilate with 
each other simultaneously, i.e. ko = ki = 0. 



FIG. 2. (color online) Localization length normalized by the 
system size, A, as a function of disorder strength for \^\ = 0.45 
and various system size L. For a guide to the eyes, we put 
three vertical dotted lines for the scale-invariant critical points 
where the localization length barely changes with L. We note 
that the system is highly anisotropic, and the transfer ma¬ 
trix calculation along x-direction (in-plane direction) seems 
to suffer large corrections to scaling. 


Localization Length and Phase Diagram. — By the TM 
method, the localizaiton length along the 2 ;-direction (the 
stacking direction) is calculated for various system sizes 
(N = L‘^ X 10^ - 10^ with L = 8,10,12,14, L being 
the linear dimension of the cross section) as a function of 
W and \/3\. In the presence of finite interlayer coupling 
(I yd I 7 ^ 0), the QH critical point in the 2D limit becomes 
a finite window of a metallic region between AI and CL 
With larger coupling, the finite region of the metallic 
phase become two distinct metallic phases, which are sep¬ 
arated by a critical line. One is in stronger disorder side 
and is connected with the 2D QH critical point, while the 
other is in weaker disorder side and belongs to the same 
phase as the 3D WSM phase in the clean limit (Fig. [1]). 
These two metallic phases are always separated by a tran¬ 
sition line where the localization length ^ normalized by 
the linear dimension of the system size L barely changes 
as a function of L (Fig.[2j) |2ll423|.The scale-invariant be¬ 
haviour of A = ^/L suggests existence of quantum critical 
line between these two metallic phases. 

Density of States. —To characterize this critical line, 
we computed the density of states (DOS) in terms of 
the KPE method for several different cubic system sizes 
{N = L^ with L = 40, 48, 60 and 80), to extrapo¬ 
late the thermodynamic behaviour of the DOS based on 
Pl{E, W) = p{E, W) + b{E, W)/L^ 0. Here pl{E, W) 
denotes the DOS as a function of the electron energy E 
and W for the linear dimension L. The intercept p{E, W) 
as a function of l/L^ can be regarded as the density of 
states in the thermodynamic limit. Fig l3] shows the DOS 
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near the zero energy for those disorder strengths around 
the phase transition between the two metallic phases. 
The DOS for the weaker disorder side vanishes at the zero 
energy, where p{E^ W) becomes a parabolic function in 
E. The DOS in the stronger disorder side acquires a finite 
value at E" = 0 This feature is consistent with 

self-consistent Born analyses (SCBA) 28 . 2^ and 

renormalization group (RG) analysis |2J, [30|, suggesting 
that the former metallic phase belongs to the same phase 
as the WSM phase in the clean limit; renormalized WSM 
phase whose zero-energy state (an electronic state at the 
band touching point; E = 0) has an infinite life time. 
Meanwhile the latter phase is characterized by the zero- 
energy state having a finite life time; diffusive metallic 
(DM) phase. 


In the renoarmlized WSM phase, the density of states 
near E = 0 is determined by an effective Hamiltonian lin¬ 
earized near the gapless points, Ek = where 

V is a renormalized velocity and p = k — ko the momen¬ 
tum distance from the respective Weyl point ko. Each 
Weyl point contributes to DOS near E = 0 as p{E) = 

^ be- 

ing an averaged velocity. On increasing the disorder 
strength, the renormalized velocity v evaluated from 
the coefficient of E^ decreases toward a certain critical 
point {W = Wc)] purple squares in Fig. [Sb- At the 
same point, p(0) starts to take a finite non-negligible 
value; blue circles in Fig. [3b, where the SCBA identi¬ 
fies p(0) as inversely proportional to the life time of the 
zero-energy state [20|. We also find that all the data 
points for the density of states near E = 0 collapse 
into a single-parameter scaling function fo rm. p{E) = 
f{\E\5-^''), with S =\W- We\/Wc IlliH. By 
fitting the DOS curves into this scaling function, we de¬ 
termine the dynamical exponent z = 1.53 ± 0.03 and 
z/ = 0.84 ± 0.1. [ 2 ^ The scaling analysis in combination 
with the SCBA unambiguously conclude the existence of 
the quantum critical line between DM phase and WSM 
phase at W = Wc- 


Two-Terminal Conductance and Conductivity — To 
characterize the renormalized WSM and DM phases, we 
calculated the two-terminal conductance G in the in¬ 
plane direction (x-direction) for the cubic system size 
with various linear dimensions (E^ with E = 6, 8, • • • , 30) 
via the transfer matrix method [32|. In the other two 
spatial directions {y and 2 ;-directions), we impose either 
periodic boundary condition (G^) or open boundary con¬ 
dition (G^). In the Chern insulator (Cl) phase, G^ is 
nearly zero at E = 30, while, for the two metallic phases, 
an in-plane bulk conductivity is obtained from G^ by 
the linear fitting (G^ = a^L -f b). 

The calculated result shows that finite conductivity in 
the DM and the WSM phase reduces to zero toward W = 
Wc (Fig. 11]). When combined with the SCBA result, this 
behaviour is consistent with the Einstein relation for the 


(a) 




FIG. 3. (color online) (a) Density of states (DOS) at |/3| = 0.6 
near the zero energy state for those W near the critical point 
intervening between DM and WSM phase. W increases from 
bottom to top. (b) Zero-energy density of states p{E = 0) and 
cube of averaged velocity F as a function of W. The latter is 
evaluated from the fitting of p(E) to a -h bE^ with ex 1/b 
and free parameter a (see also text). The vertical dot-dashed 
line and dashed line denote the critical disorder strength at 
which V and p(E = 0) vanish, respectively (see Ref. for 
their determinations). 


conductivity cr(E); 

aiE) = e^D{E)p{E). (2) 

The diffusion constant D{E) is given by the life time r 
and the averaged velocity v as D{E) = SCBA re¬ 
lates p{E) with the inverse of the life time as p{E) = 
[ 2 ^. Thus the relation tells that cr(E = 0) = 
Sc^v‘^I{'kW‘^) remains finite not only in the DM phase but 
also in the renormalized WSM phase. Moreover, since v 
vanishes at the critical point as above [ 2 ^, the Einstein 
relation also dictates that the conductivity reduces to 
zero at the critical point, being consistent with our nu¬ 
merical observation (blue points near IT = 2 in Fig. |4|). 
Note also that cr(E = 0) vanishes as from the DM 

side while v vanishes as from the WSM side ( 2 ^ . 
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FIG. 4. (color online) Bulk conductivity and conductances of 
the cubic system (L^) at |/3| = 0.45 (inset: at \I3\ = 0.50) as 
a function of disorder strength W\ the bulk conductivity cr^ 
(blue points with line), conductance with the periodic bound¬ 
ary condition (L = 30; red points with line) and that with 
open boundary condition (L = 30; yellow points with 
line). The vertical dotted lines around W — 1.92 (|/3| = 0.45) 
and W = 1.96 (|/3| = 0.50) stand for the WSM-DM transi¬ 
tion points determined by the DOS, while the vertical dotted 
lines around W = 8.7 (|/3| = 0.45) and W = 8.9 {\^\ = 0.50) 
and that around W = 1.35 (|/3| = 0.45) are the DM-AI tran¬ 
sition points and CI-WSM transition point determined from 
the localization length, respectively. 


With our expression for the conductivity a{E = 0) oc 
we see that the conductivity vanishes continuously from 
the WSM side as a{E = 0 ) oc This is qualita¬ 

tively consistent with Ref. [H, which predicts a contin¬ 
uously vanishing conductivity with a different exponent 
cr(F^ = 0) oc (5. If we assume that z = (i/ 2 , our exponent 
in the WSM side coincides with that from the DM side. 

and signficantly differ from each other in the 
Chern insulator phase and the renormalized WSM phase 
(Fig. |4]). The difference can be attributed to the surface 
conductance due to the chiral surface states. In the renor¬ 
malized WSM phase, the zero-energy state has an infinite 
life time r, so that the SCBA is fully justified; neglected 
Feynman diagrams in the approximation are smaller than 
those included at least by a factor 1/r. The SCBA re¬ 
sult shows that the Green’s function after quenched im¬ 
purity averaging has either three or two pairs of Weyl 
points in the momentum space (inset of Fig. [T]). In the 
region-II WSM phase, two pairs of Weyl points appear 
qX k = ( 0 , 7 r, 7 r ± ko) and ( 7 r, 0 , 7 r ± /cq), while another 
one pair appears at (7r,7r, =F^i) the region-III WSM 
phase. These Weyl points in the bulk electronic state re¬ 
sult in a surface chiral Fermi arc state with left-moving 
chirality connecting from the surface crystal momentum 
^z) = (O 5 tt — /co) to ( 0 , —TT + /co), and two other chi¬ 
ral arc states with right-moving chirality, one connecting 


FIG. 5. (color online) AG = G° — G^ for |/3| = 0.45 (inset 
for 1^1 = 0.55), where G° and G^ are calculated for the cubic 
system (L^ = 30^). The red solid line is a surface conductance 
evaluated from eq. along with the SGBA result for the 
same parameter set. 


from (tt, tt — k^) to (tt, ki) and the other from (tt, —ki) to 
(tt, —tt E ko). They in total lead to a finite two-terminal 
surface conductance which depends on the length of each 
Fermi arc. The left-moving chiral arc and right-moving 
arc cancel each other by the intra-surface backward scat¬ 
tering due to impurities, so that the surface conductance 
is expected to be expressed with the difference between 
the length of right-moving arc (27r — 2ko — 2ki) and that 
of the left-moving arc {2ko); 


‘2‘k{} k\)L 

h TT 


( 3 ) 


In fact, ko and ki obtained from the SCBA along with 
eq. ([3j) reproduces AG from the numerics at the quan¬ 
titative level for W < Wc (Fig. [5]). Importantly, paired 
two Weyl nodes do not annihilate with each other even 
at VF = Wc] the associated chiral surface Fermi arc state 
does survive up to VF = IFc. For W > VFc, each Weyl 
node acquires a finite line broadening (inverse of life¬ 
time), so that the chiral Fermi arc states near the two 
ends of the arc start to be mixed with bulk states and 
lose their protected surface conduction property. Mean¬ 
while, those arc states away from the two ends still con¬ 
tribute to robust surface conductions for W > Wc- 0 
This picture is consistent with our numerical observation 
of the finite surface conductance AG for W > IFc, which 
is smaller than the SCBA estimate from Eq. m in the 
DM phase. 

Gonclusion — In this Letter, we introduce a simple lat¬ 
tice model in the unitary class, which enables us to study 
WSM phase, disorder-induced DM phase and quantum 
phase transition between these two phases. The criti¬ 
cal nature of the phase transition is confirmed by the 
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scale-invariant behaviour of the localization length and 
by the scaling analysis of the density of states. The calcu¬ 
lated critical exponents take those values close to critical 
exponents previously obtained in a lattice model of the 
symplectic class [26|. We show that longitudinal con¬ 
ductivity at the zero-energy state remains finite in the 
renormalized WSM phase as well as in the DM phase, 
while it reduces to zero toward the transition point. With 
the help of the SCBA, we show that this observation is 
consistent with the Einstein relation for the conductivity 
and Wegner’s conductivity scaling [s^. The two-terminal 
conductance shows a significant residual surface conduc¬ 
tance at the transition point, which can be attributed to 
the “renormalized” chiral Fermi arc states. 

Acknowledge The authors would like to thank Dr. Koji 
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ported by JSPS KAKENHI Grants No.l5H03700 and 
No. 24000013 and by NBRP of China (Grant No. 
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Note added. - Recently we became aware of indepen¬ 
dent numerical works [HI that treat similar situation but 
focus on a different aspect. 

SUPPLEMENTAL MATERIALS 

Weyl semimetal phases and chiral Fermi arc in 
layered Chern insulator 

The tight-binding Hamiltonian for layered Chern insu¬ 
lator [3, Q reduces to the following 2x2 Hamiltonian in 
the momentum space, 

H(k) = aoCTo + a • (T (4) 

with a = {cFx^cFy^crz) are Pauli matrices and 

a»(fc) = + (tp - f.Xcost. +co.t.) - ((', + 

^ - {tp{cos kx ^ cosky) - {t'g -tp)cosk, 

ai{k) = —2tspSmky, 
a2{k) = -2tspAiikx, 

respectively. For simplicity, we take tp = tg > 0 and 
t's = —t'p so that the /c-dependence of uq can be omit¬ 
ted and the direct band gap always becomes the global 
band gap A(/c) = 2 min|a(/c)|. If any, the gap clos¬ 
ing occurs at the high symmetric k lines, i.e. k = 

( 0 , 0 , kz)^ (tt, 0 , kz)^ ( 0 , TT, kz)^ (tt, tt, kz). At these points, 
the gap is given as a function of kz 



FIG. 6. (color online) Phase diagram of layered Chern in¬ 
sulator in the clean limit. 01 and Cl denote ordinary band 
insulator and Chern insulator respectively, while WSM de¬ 
notes the Weyl semimetallic phases. The explanation for four 
distinct regions of the WSM phase is given in the text of this 
supplemental material. 

in Fig. [6l The phase diagram comprises of Chern insula¬ 
tor phases and conventional band insulator phases, and 
four Weyl semimetal phases; WSM in the region I has 
one pair of two Wely points having opposite magnetic 
charges (either dX k = (0,0, ±*) or (7r,7r, ±*)), WSM in 
the regions H has two pairs (at k = (tt, 0, ±*), (0, tt, ±*)). 
In the region HI and VI, it has three pairs (either at 
k = (0,0, ±*), (7r,0,±*), (0,7r, ±*), or at (7r,7r, ±*), 
(7r,0,±*), (0,7r, ±*)) and four pairs respectively. In the 
Inain text, we took a = — so that, on increasing |/3|, 
the system enters the region HI first (^ < |yd| < |) and 
’then region IV (| < |/3|). Figure [3 shows surface chi¬ 
ral Fermi arcs in the case of |/3| = 0.6, where the right- 
moving chiral Fermi arc connects from a surface momen¬ 
tum point {kx^kz) = (0,7r — k^) to (0, —tt -h k^) with 
ko = cos“^[^], while the two left-moving chiral Fermi 
arcs connect from (kx^kz) = (7r,/ci) to {tt^tt — ko) and 
from {kx, kz) = (tt, —ki) to (tt, -ir-^-ko) respectively with 
ki = cos“^[— 


Self-consistent Born analysis 


A{k) = 2{ts + tp) 


a ■ 


cos kx + cos ky 


ti-tt 


-h P cos kz 


a = and /3 = = 2 *. • 


This 


A single-point Green function is averaged over the on¬ 
site impurity potentials. 


zts \r^,(h O = IIh n/\ 

gives a phase diagram subtended by a and (3 depicted )\a,/3 — 


1 


E-Ho-V±iS 


^ „p(5) 














6 



FIG. 7. (color online) (a) Energy dispersion of chiral Fermi 
surface states localized at one surface (red) and at the other 
(blue). Those energy states with grey color points are bulk 
states, (b) Location of chiral Fermi arcs are depicted on the 
surface crystal momentum space (the kx-kz plane), where a 
red solid line stands for a right-moving chiral Fermi surface 
state [dEjdkx > 0), and red dotted lines are for left-moving 
chiral Fermi surface states {dEjdkx < 0). 


where \k,a) = \k)\a), \k) a plane-wave state with mo¬ 
mentum /c, and \a) either s or p orbital. ELo is the non¬ 
interacting Hamiltonian in the clean limit and V denotes 
the impurity potential part in Eq. ( 1 ) in the main text, 
i.e. V = + Ex^'p(®) 4 ,pCa!,p- The im- 


purity averages are taken as follows, 

(^Vs{Xj)Vs{Xm))^^^ = (^Vp{Xj)Vp{XTji)'^^^^ 

\ _ 

dv 1 = Sj^rn ^2 ~ 

with {vp{xj)vs{xm))imp = 0. Within the self-consistent 
Born approximation, the Green function after quenched 
average takes a form of [G±{k, k')]a^p = [G±{k)]a^p 

where the 2 by 2 matrix G±{k) is given by itself, 

Gl\k) = GoXik) -YJfEE <^f^G±{q)cr^ 

q /Lt=0,3 

with N the number of sites. GQ^{k) denotes a non¬ 
interacting Green function, Gq ^(/c) = {uj ± iS) — H{k) 
with Eq. dH. The approximation takes into account 
all the ‘non-crossing’ diagrams, while ‘crossing’ diagrams 
(Feynman diagrams containing crossings of two impurity 
potential lines) are smaller than non-crossing diagrams 
by an additional factor proportional to . Thus, the 
approximation is justified in the limit of long life-time 
quasiparticle (r“^ ^ 0 ), which is indeed the case with 
the zero-energy state {E = 0) in the renormalized Weyl 
semimetal (WSM) phase (see below). 

The averaged Green function is expanded in terms of 
the Pauli matrices with complex-valued coefficients, 

3 

Gl^{k) = bocro (6) 

j = l 

The symmetric treatment of the impurity average ( 
{vsVs)imp = {vpVp)imp and {vsVp)imp = 0 ) cuablcs hi{k) = 
ai{k) and 62 (^) = Ci 2 (k). bp = E iS — jo and 
bs{k) = as{k) +73 with A^-independent complex-valued 
coefficients 70 and 73 . These two are determined by the 
following coupled self-consistent equations; 

_ ^ _ 7o ~ _ 

^ k 4(^) + 4(fe) + (a3(fe) + 73 )^ - {E+ - 70 )^ ’ 

( 7 ) 

_ Tf a3(fe)+73 

alik) + ai(fe) + (a3(fc) + 73 )^ - iE+ - 70 ) 2 ' 

(8) 



with E_p = E iS 

For the zero-energy state {E = 0 ), Eqs. ([ 7 )) and m 
reduce to a simpler set of coupled equations. The simpler 
set of equations are only for Imyo = and Re73 = M3 
while Reqo = 0 and Im73 = 0; 


1 

r 


1 K 


'n 




1 

—2 . . ’ 

( 9 ) 

E (fe) +T-2 


03 (fe) + M3 

_2 ’ 

(10) 

E {k) + T-^ 


al{k) + (a3(fc) + M3)2. 

(11) 





















7 


For weaker disorder case, the first equation can be satis¬ 
fied only by putting =0. The solution with = 0 
(zero-energy state with infinite life time) corresponds to 
either the renormalized Weyl semimetal (WSM) phase or 
gapped Chern insulator (Cl) phase. A critical disorder 
strength Kc above which the coupled equations can have 
a solution with finite life time^“^ ^ 0 ) is determined 
by the following gap equation [25|, l28| , 


i = —E 

N 




( 12 ) 


where M 3 in E{k) in the right hand side is obtained from 
Eq. ([To]) with =0. In fact, for K > Kc, finite r can 
satisfy 


1 


K y. 1 
N\^E\k)+T-^' 


(13) 


The finite-r solution (zero-energy state with finite 
life time) corresponds to the diffusive metallic (DM) 


phase [25 


For K < Kc, we solve Eq. (HHD with r ^ = 0 in favor 
for M 3 (renormalized WSM or Cl phases). For K > Kc, 
we solve Eqs. ([TQ|) and ([T3|) for M 3 and r (DM phase). In 
either cases, an electronic dispersion is renormalized by 
the disorder-induced mass term M 3 . In the present set 
of parameters with a = —^, M 3 becomes positive. Fig. [ 6 ] 
suggests that the positive mass term with a = —^ drives 
the Cl phase and the region-III WSM phase (with three 
pairs of Weyl points) near |/3| = 0.5 into the region-II 
WSM phase (with two pairs of Weyl points). Namely, 
for the Cl phase {P > —^ = a), larger positive mass M 3 
realizes a -h M 3 > (3, where two pairs of Weyl points are 
created ni k = (tt, 0 , tt ± ko) and ( 0 , tt, tt ± ko) with 


ko = cos [- - -J 


(14) 


For the region-III WSM phase {^ < —\ = —1 — a), n 
positive mass term can satisfy 1 -h cr + M 3 > —/3, so that 
a pair of two Weyl points nt k = ( 7 r, 7 r, ±/ci) annihilate 
with each other with 


ki 


= cos 


l + a +Ms 

P 


(15) 


into Eq. (3) in the main text, we calculate the surface 
conductance due to the “renormalized” chiral Fermi arc, 
AG = Go — Gp, as a function of the disorder strength 
(Fig. 5 in the main text). For K < Kc, AG thus ob¬ 
tained quantitatively reproduces the numerics. 

Importantly, the SCB analysis shows that paired two 
Weyl nodes do not annihilate with each other at the 
semimetal-metal quantum phase transition point {K = 
Kc). Thereby, the associated chiral surface Fermi arc 
state also does survive at the transition point. Instead 
of annihilating with each other, the two nodes acquire 
a finite line broadening individually', an inverse of their 
life time (t“^) evolves continuously from zero to a finite 
value for K > Kc (as described above). Due to this line 
broadening of the Weyl nodes, those chiral surface states 
near the two ends of the arc start being mixed with bulk 
states, losing their protected surface conduction nature. 
On the one hand, those chiral surface states far away 
from the two ends still remain intact for K > Kc, be¬ 
cause they are well separated from bulk states (Fig. [ 8 |). 
The line broadening of the Weyl nodes evolve continu¬ 
ously, so that the latter surface states can still contribute 
to a residual surface conduction AG for those K moder¬ 
ately greater than Kc. In fact, the numerics (Fig. 5 in 
the main text) shows that a surface conductance remains 
finite for W < 2Wc, while the DM phase itself ranges up 
to a much larger disorder strength. 

Note also that the persistence of paired two Weyl nodes 
ni K = Kc is consistent with the vanishing effective ve¬ 
locity V at the transition point (observed by the density 
of states analysis). The SCB analysis explains T = 0 
nt K = Kc by an anomalous scaling of a renormalized 
chemical potential. Specifically, the SCB analysis shows 
that, for smaller E*, the real part of 60 in Eq. m is scaled 
with E for K < Kc and with \/E for K = Kc. 

Within the self-consistent Born approximation, the 
density of states is proportional to an inverse of the life 
time r. 


1 1 


1 




E-Ho-V + iS 
2 


tttK 


(16) 


By solving Eq. m, we determine a phase boundary be¬ 
tween diffusive metallic phase and renormalized WSM 
phases. By solving Eq. m with r ^ = 0 in favor 
for Ms, we determine a boundary between the region- 

II WSM phase and Chern insulator from a + M 3 = /3, 
and a boundary between the region-II WSM and region- 

III WSM phases from l-\- a-\- Ms = —p. This completes 
a phase diagram obtained from the self-consistent Born 
analysis (inset of Fig. 1 in the main text). By solving 
Eqs. (|9|) and (flO]) for M 3 and r, we also obtain Eqs. m 
and m-, “renormalized” locations of the Weyl points 
in the momentum space. By substituting ko and ki 


Numerical solutions of Eqs. © and dHj) dictate that r ^ 
for smaller E is proportional to E^ in the renormalized 
WSM phases {K < Kc) [25|,[28|, while being scaled to \/E 
at the quantum critical point {K = Kc). The solutions 
also find that the inverse of the life time of the zero- 
energy state in the DM phase {K > Kc) grows linearly in 


K — Kc 23 , l28| . r ^ (X for E < Ec is consistent with 


the behaviour of p{E) for W < Wc obtained from the 
kernel polynomial expansion (Fig. 3 in main text). oc 
y/E at K = Kc results in a dynamical exponent 2)scb = 2, 
which is different from the exponent evaluated from the 
numerics {z 1.5). (x K — Kc at E = 0 for K > Kc 
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W>Wc (K>Kc) 

FIG. 8. (color online) Schematic pictures of an energy- 
momentum (surface crystal momenta) dispersion for bulk 
states with a pair of two Weyl nodes (black) and for corre¬ 
sponding chiral surface Fermi arc state at = 0 (a red line), 
(a) K <Kc and (b) K > Kc. As in Fig. [3 the surface crystal 
momenta comprise of two momenta, kx and kz^ while kx is 
omitted for simplicity; the surface Fermi arc state has a chiral 
dispersion along the kx direction. Note also that we only show 
one pair of two Weyl nodes for clarity. The disorder endows 
bulk states with a finite life time (line broadening), whose 
effect is represented by a green hatch in the figure. Within 
the green hatch, the crystal momentum is not well-defined 
even after quenched averaging of disorder potentials; any two 
states within the hatch are mixed with each other and indis¬ 
tinguishable due to the disorder, (a) for K < Kc, sm inverse 
of the life time (line broadening) for the two Weyl nodes is 
zero, so that all the chiral arc states at = 0 connecting 
these two nodes provide protected surface conductions, (b) 
for K > Kc, those chiral fermi arc states at F^ = 0 which 
are close to the two ends of the arc (those red dotted lines 
covered by the green hatch) are mixed with bulk states due 
to the finite line broadening of the two Weyl nodes; they will 
not contribute to surface conductions. Meanwhile, those arc 
states at F^ = 0 which are well separated from the two ends 
(red solid line not covered by the green hatch) are not mixed 
with bulk states, providing protected surface conductions. 


in combination with the dynamical exponent Zgcb = 2 
gives the other exponent z^scb = 1- Discrepancies of the 
critical exponents obtained from these two calculations 
stem from the neglected Feynman diagrams (‘crossing’ 
diagrams) in the self-consistent Born approximation. 

The Einstein’s relation relates the longitudinal electric 
conductivity cf{E) with the diffusion constant D{E) and 
the density of states p{E), 

a{E) = e‘^D{E)p{E). (17) 

The diffusion constant D{E) is given by mean-free length 
I = vr and life time (mean free time) r as D{E) = 
P/3r = v^Tj3, where v stands for an effective velocity 
of electrons. As shown above, the density of states p{E) 
is inversely proportional to the life time within the self- 
consistent Born approximation. Thus, when the energy 
is set to the zero in the renormalized WSM phase, p{E^ 
vanishes as 1/r, D{E) diverges as r, and the conductivity 
remains constant; a{E = 0) = 2v‘^/{S ttK). This dictates 
that the conductivity at E = 0 remains finite not only 
in the diffusive metallic phase but also in the renormal¬ 
ized WSM phase, while vanishing at the quantum critical 
point intervening these two phases where T ^ 0 (see the 
main text). 


Chiral surface states in layered Chern insulator and 
in renormalized Weyl semimetal phase 

Eigenstates of Eq. (1) in the main text for a given 
disorder realization are numerically calculated with the 
periodic boundary condition along the chiral (x) and the 
stacking {z) directions, and open boundary condition for 
the other directions (y). System size L = SO, and numer¬ 
ical calculation has been done by sparse matrix diagonal- 
ization algorithm of Intel MKL/EEAST. Eigure [9] shows 
spatial density distributions of an eigenstate u, 

Pu{x,y,z) = \^pu,six,y,z)f + \tpt,,p{x,y, zp . 

We took the states u which are closest to the zero en¬ 
ergy {E = 0) in the layered Chern insulator phase region 
(a) {\/3\ = 0.45 and W = 0.8), and in the renormal¬ 
ized Weyl semimetal phase region (b) {\/3\ = 0.45 and 
W = 1.7). In both cases, the eigenstates are nearly lo¬ 
calized at the boundaries {y = ELj2), suggesting that 
zero-energy state or eigenstates close to the zero energy 
are usually chiral surface states. Eigure [9] also indicates 
that the chiral surface state in the renormalized WSM 
phase is spatially extended within surfaces, while the chi¬ 
ral state in the layered Cl phase is extended only along 
the chiral direction and localized along the stacking di¬ 
rection. 
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(a) 


(b) 


FIG. 9. (color online) Wave function density examples for 
layered Chern insulator (a) (/3 = 0.45, VF = 0.8, = 

—0.00000859 •••), and a wave function density for Weyl 
semimetal (b) (/3 = 0.45, W = 1.7, = -0.0000688 • • • )• In 

the former case, the charge density appears one dimensionally 
along X direction on the x — z plane, while in the latter case, 
the charge density appears two dimensionally on the x — z 
plane. 
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Finite size scaling of the density of states 

The density of states (DOS) is calculated for differ¬ 
ent system sizes (L = 40, 48, 60, 80) in terms of the 
kernel polynomial expansion method [l9|. The calcu¬ 
lated DOS is fitted with the finite size scaling form, 
Pl(B,W) = a(E,W) + b{E,W)/L^ which gives 
the DOS in the thermodynamic limit as the intercept 
a{E, W). An example of the fitting is shown in Fig. [10] 
(for P = —0.6, E = 0.0 and W = 1.5 < Wc). The error 
bar obtained from this finite-size-scaling fit is shown in 
Fig. [TT] for typical values of E and W with (3 = —0.6. 
The error is estimated with 95% confidence (2cr). Since 
the fitting has only two degrees of freedom, the estimated 
error bar is relatively large. 


Parabolic fitting and the velocity 


FIG. 10. (color online) Finite size scaling of the DOS at = 0 
and VF = 1.5. The blue squares denote numerical data of the 
DOS for four different system sizes, and the red dashed line 
is the linear fit to 1/L^. 



-0.15 -0.1 -0.05 0 0.05 0.1 

E 

FIG. 11. (color online) The DOS from the finite size scaling 
for typical values of E and IF. The 95% confidence intervals 
are indicated as error bars . 


In the renormalized WSM phase (IF < IFc), the av¬ 
eraged velocity of the Weyl point is obtained from a fit¬ 
ting of the DOS near F^ = 0 by a parabolic function 
p{E) = a -f hE‘^. Fig. [12] shows an example of the fit¬ 
ting {/3 = —0.6, IF = 1.5 < IFc). On increasing \E\, 
the DOS curve changes from the F^-square behavior to a 
F^-linear behaviour. When IF becomes closer to IFc from 
below, the F^-square region becomes narrower in energy, 
indicating the reduction of the velocity. To estimate the 
curvature only around F^ = 0 from the fitting, we first 
observe this crossover by eye and then introduce a cut¬ 
off for p{E) to exclude those data points in the E'-linear 
region (For example, p{E) < 0.013 in Fig. [T2|). 


Determination of the critical point 

The phase boundary between WSM and DM is deter¬ 
mined by an onset of finite DOS at F^ = 0 and also by an 
onset of finite averaged velocity of the Weyl cone. The 
DOS at F^ = 0 vanishes for smaller disorder strengths. It 
starts to take a finite value when IF is greater than a crit¬ 
ical disorder strength. To determine this critical disorder 
strength without ambiguity, we first select a set of data 
points in the p(0) — IF figure (typically five to eight data 
points), which are in the DM region within the finite- 
size-scaling error bar (explained above), but which are 
sufficiently close to the WSM region. We applied a IF- 
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FIG. 12. (color online) A parabolic fitting of the DOS at 
W = 1.5 and [3 = —0.6. The blue solid line is from numerical 
data and the red dashed line is a F^^-fitting curve. To excludes 
those data points away from = 0, we introduce a cutoff for 
the DOS {p < 0.013), which is depicted by a black dashed 
horizontal line. 


linear fitting for these data points, to regard the intercept 
of the fitting curve with the /)(0) = 0 axis as the critical 
disorder strength (Fig. [T^a)). 

On increasing the disorder strength in WSM phase, the 
inverse of the parabolic coefficient 1/6 decreases continu¬ 
ously, only to converge to a small but non-zero constant 
value Cres in dm phase side. The small non-vanishing 
Cres is due to non-zero error bar stemming from the finite 
size scaling and the finite truncation of KPE. To derive a 
critical disorder strength as an onset disorder strength of 
finite averaged velocity, we regard this residual value C^es 
to be effectively zero and determine the critical strength 
from an intersect of a curve for 1/6 as a function of W 
with 1/6 = Cres- To be more specific, we firstly fit out 
Cres from a set of those data points in the DM region 
which are close to the WSM region (data points encom¬ 
passed by a grey dotted circle in Fig. [T3l b)). We then 
collects another group of data points from the WSM re¬ 
gion which are sufficiently close to the DM region (data 
points encompassed by an orange-colored dotted circle 
in Fig. [T3f b)). As above, we applied a VF-linear fitting 
for the latter group of the data points, and regard its in¬ 
tercept with 1/6 = Cres as the critical disorder strength 
(Fig.imb)). 


Single parameter scaling 

According to a scaling argument [ 2 ^ [HI , the density of 
states near the zero energy state follows scaling functions 
near the quantum critical point; 

p{E) = (18) 


(a) 



(b) 
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FIG. 13. (color online) (a) Determination of the critical dis¬ 
order from the density of states (DOS) at E = 0. We 

selected a group of data points in the diffusive metal (DM) 
side and tried a linear fitting (red dashed line). The inter¬ 
cept with the p(0) = 0 axis is identified with the critical 
disorder strength . (b) Determination of the critical 

disorder from the averaged velocity of the Weyl cone 

1/6. The data points in the DM side, which are near the 
WSM phase, take approximately small but constant value 
Cres (black dashed line). We applied a linear fitting to a 
group of data points in the WSM side (red dashed line). The 
intercept with 1/6 = Gres is regarded as the critical disorder 
strength , below which the averaged velocity starts to 

take a finite value. 


where a characteristic length ^ scaled with the distance 
measured from the critical disorder strength, 6 = \W — 
Wc\/Wc with the critical exponent z/ as ^ oc . z de¬ 
notes the dynamical exponent and the spatial dimen¬ 
sion d is 3 in the present case. The critical disorder 
strength Wc is chosen to be the mean value between 
W^SM ^ since these two coincides with each 

other within the error bar. A scaling function /(•••) for 
the WSM phase and that for the DM phase are generally 
different from each other, while the critical exponents 
(z/ and z) are universal. At the critical point, the 5- 
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dependences in Eq. m should cancel exactly, leading to 
piE) oc = \E\^d-z)/z_ 

lated DOS at VD = Wc is nearly linear in \E\, suggesting 
that 0 ~ d/2 SSI. A comparison between the numer¬ 
ical data and the scaling form gives 2 ; = 1.53 ± 0.03. 
In WSM phase, the DOS near = 0 is scaled with 
\E\^~^^ so that the scaling form requires that its coef¬ 
ficient should be scaled with p{E) oc 
Fitting of an coefficient of E^ from Fig. 3b in the main 
text along this scaling form, i.e. 1/5 oc ^ gives 

out d{z — l)iy = 1.24 ± 0.06. In DM phase, the DOS at 
= 0 is scaled with p(0) = along which 

the DOS at F = 0 from Fig. 3b in the main text is fit¬ 
ted, giving {d — z)u = 1.47 ± 0.06. These two fittings in 
combination with 2 ; = 1.53 ± 0.03 give us 

^WSM ^ 0.78 ±0.08, (19) 

^DM ^ 1.00 ±0.13, (20) 

respectively. Though they do not coincide exactly, the 
discrepancy is at the same order of our estimated error 
bar, giving a weighted average u as 

z/= 0.84 ±0.10. (21) 

In fact, with this critical exponent and 2 ; = 1.53 ± 0.03 
in hand, we can successfully fit all the data points of our 
DOS for p{E) > pc into two scaling functions, provided 
that pc is chosen on the same order of error bars stem¬ 
ming from the finite size scaling (see above). As shown in 
a log-log plot of versus \E\6^^ (Fig. [14]), all the 

data points collapse into two branches, where one branch 
(lower branch) corresponds to a scaling function for the 
WSM phase and the other (upper branch) to a scaling 
function for the DM phase respectively. With these re¬ 
sults, we conclude that the critical nature of the WSM- 
DM phase transition in the unitary class follows the scal¬ 
ing law of Eq. (p±|) with 2 ; = 1.53±0.03 and u = 0.84±0.1. 



ln(|E|6-"") 


FIG. 14. (color online) A single parameter scaling of the 
density of states for \/3\ = 0.6 and W near the critical point 
with p{E) > 0.02. All the data points collapse onto two 
distinct branches which correspond to a scaling function for 
WSM phase (lower branch) and that for DM phase (upper 
branch), respectively. 


a function of the cubic system size L, while the conduc¬ 
tance in the WSM phase for \f3\ = 0.45 and 0.50 show 
much less prominent oscillation. The larger oscillations 
for I yd I = 0.55 and 0.6 make it hard to determine the 
bulk conductivity in the WSM phase with smaller error 
bar (not shown in Fig. [T5|). The oscillation is due to (i) 
the very long mean-free path in the WSM phase and/or 
due to (ii) a mismatch between (A^ arbitrary integer 
less than L) and ko (or ki) (the renormalized positions 
of the Weyl cones along the kz axis). 


Conductance and conductivity at |/3| = 0.55,0.6 

The two-terminal conductance and conductivity for 
I yd I =0.55 and 0.6 are shown in Fig. [151 The conductance 
along the x direction (in-plane direction) are calculated 
for various system size {L = 6, 8,10, • • • , 30), with peri¬ 
odic (G^) or open boundary condition (G^) in the y and 
2 ; directions. We prepared 20 uncorrelated disorder real¬ 
izations, over which the calculated conductance are aver¬ 
aged. The in-plane bulk conductivity is obtained from 
G^ by the linear fitting (G^ = a^L ± c). We found that 
a finite bulk conductivity in the DM phase continuously 
reduces to zero toward the WSM-DM transition point. 
This behaviour is consistent with the self-consistent Born 
analyses 0 and a generalized Wegner’s scaling relation 
[ 21 , 0 . 

Note also that the two-terminal conductance in the 
WSM region for |yd| = 0.55 and 0.6 strongly oscillate as 
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